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ABSTRACT 

We present a survey of the X-ray emitting ejecta in the Cassiopeia A su- 
pernova remnant based on an extensive analysis of over 6000 spectral regions 
extracted on 2.5-10" angular scales using the Chandra 1 Ms observation. We 
interpret these results in the context of hydro dynamical models for the evolution 
of the remnant. The distributions of fitted temperature and ionization age are 
highly peaked and suggest that the ejecta were subjected to multiple secondary 
shocks. Based on the fitted emission measure and element abundances, and an 
estimate of the emitting volume, we derive masses for the X-ray emitting ejecta 
as well as showing the distribution of the mass of various elements over the rem- 
nant. The total shocked Fe mass appears to be roughly 0.14 M©, which accounts 
for nearly all of the mass expected in Fe ejecta. We find two populations of Fe 
ejecta, that associated with normal Si-burning and that associated with ct-rich 
freeze-out, with a mass ratio of approximately 2:1. Surprisingly, essentially all 
of this Fe (both components) is well outside the central regions of the SNR, pre- 
sumably having been ejected by hydrodynamic instabilities during the explosion. 

We discuss this, and its implications for the neutron star kick. 

Subject headings: hydrodynamics - ISM: individual (Cassiopeia A) - supernova 
remnants - X-rays: ISM 


1 Goddard Space Flight Center and Johns Hopkins University 

Una. Hwang- lOgsfc .nasa.gov 

2 Code 7674L, Naval Research Laboratory, Washington DC 20375 


laming@nrl . navy . mil 



- 2 - 


1. Introduction 

Theory and observations are merging toward consensus that core-collapse supernova 
explosions are intrinsically asymmetric, even if the progenitor had an initially symmetric 
configuration. Recent three-dimensional simulations show that strong dynamical interac- 
tions between the various burning shells can lead to large asymmetries in the progenitors, 
independent of the symmetry of the explosion (Arnett & Meakin 2011). Nor is the explosion 
itself likely to be spherical, with convective instabilities (e.g. Herant et al. 1992; Herant 1995) 
and instabilities at the accretion shock (e.g. Blondin et al. 2003; Foglizzo et al. 2007) all pos- 
sibly acting to create asymmetries. It is further anticipated that such asymmetries may hold 
the key to the core-collapse explosion mechanism (e.g. Ott et al. 2008; Brandt et al. 2011). 
They will be manifested not only in the distributions of the ejecta mass and velocity, but also 
in the recoil of the nascent neutron star, which is observed with measured space velocities 
of up to 1000-1500 km/s (Arzoumanian et al. 2002; Hobbs et al. 2005; Faucher-Giguere & 
Kaspi 2006). Simulations show that a neutron star kick that has a purely hydrodynamic 
origin can balance the total momentum of the anisotropic component, and is expected in 
the direction opposite the fastest shock expansion (Wongwathanarat et al. 2010; Scheck et 
al. 2004). 

Observations of core-collapse supernovae indicate a significant degree of explosion asym- 
metry via measurements such as time-dependent polarization changes or Doppler velocities 
(Mazzali et al. 2006; Leonard et al. 2006), but explosion asymmetries are also revealed in 
sufficiently young supernova remnants (SNRs) through their shock-heated ejecta. The most 
detailed view so far available of the ejecta in a young core-collapse SNR is provided by Cas- 
siopeia A, which at 330 years is the penultimate of the known Galactic supernova remnants. 
Its reverse shock has already progressed deeply into the explosively produced nucleosynthesis 
products, aided by the strong presupernova mass loss incurred by its progenitor through the 
likely action of a binary companion (Young et al. 2006). 

That Cas A was produced by an asymmetric explosion is by now well-established, and 
most recently by light echo observations which sample the explosion hundreds of years after 
the fact in different directions. These show variations in the ejecta velocities of ~ 4000 
km/s (Rest et al. 2010). Chandra X-ray observations have also played an important role 
in revealing the complexity of Cas A’s X-ray emitting ejecta. The X-ray emitting Si ejecta 
show a bipolar structure with jet-like features (Hwang et al. 2004; Vink et al. 2004; Laming 
et al. 2006) similar to that seen in optical (Fesen 2001, and references therein) and infrared 
emission (Hines et al. 2004). The X-ray ejecta spectra reveal the imprint of roughly a factor 
of two asymmetries in the deposition of explosion energy around the remnant (Laming & 
Hwang 2003), similar to the distribution of kinetic energies and positions of fast moving 
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optical knots (Hammell & Fesen 2008). Early Chandra observations identified regions domi- 
nated by emission from Fe ejecta well outside the projected location of Si emission that were 
proposed to be sites of overturn between the nucleosynthesis layers (Hughes et al. 2000). 
This conclusion is apparently supported by the X-ray Doppler measurements showing higher 
velocities of Fe compared to Si (Willingale et al. 2002). More recently, the dynamical struc- 
ture of the ejecta has been inferred in remarkable detail by the multi- wavelength (X-ray, 
optical, infrared) Doppler-shift maps compiled by DeLaney et al. (2010). These show the 
inner ejecta to be unshocked and in a flattened distribution ([Si II]), outlying matter ([Ar 
II], [Ne II]) to be arranged in rings on a spherical surface, with outflows emerging out of the 
surface of the remnant (Fe ejecta, jet, and counterjet). Their detailed analysis shows that 
the outflows of X-ray emitting Fe ejecta are not material mixed outwards from overturning 
of the ejecta layers, but rather, material that has pushed through the overlying ejecta and 
are now encircled by rings of material from outlying nucleosynthesis layers. 

Large-scale studies of the X-ray emission in Cas A have already been undertaken by 
Stage et al. (2006) and Helder & Vink (2008), but these focus on the properties of the 
nonthermal emission. The most comprehensive accounting of the X-ray emitting ejecta 
in Cas A to date was carried out by Willingale et al. (2003), who used XMM-Newton 
CCD observations to examine 225 regions of fixed 20" size across Cas A (Willingale et 
al. 2002). They fitted a two-component thermal model throughout to represent hot shocked 
circumstellar material and cooler shocked ejecta and inferred a total ejecta mass of 2.2 M 0 . 
They interpret the Fe K emission as forming a bipolar double cone that is associated with 
ejecta bullets that have broken beyond the forward shock. 

A detailed accounting of the X-ray emitting ejecta on fine angular scales was one of the 
primary goals of a 1 Ms Chandra observation of Cassiopeia A that was obtained in 2004, 
and is the subject of this paper. Our spectral survey includes over 6000 spectral regions 
extracted on 2.5-10" angular scales. It is distinguished from previous Chandra studies of the 
ejecta in Cas A, which have relied either on spectral imaging or the detailed spectral analysis 
of a relatively small number of regions. We are able to cover a larger area of the remnant 
with much higher sensitivity and angular resolution than was available to Willingale et al. 
(2002). In particular, we are able to show the distribution of the X-ray emitting Fe ejecta 
in detail. The Fe ejecta, being produced just outside the collapsing core of the supernova, 
provide valuable clues to understand important details of the explosion, including the recoil 
of the ejecta with the neutron star. 

We interpret the spectra in the context of the hydrodynamical models applied by us 
previously to Cas A (Laming & Hwang 2003; Hwang & Laming 2003; Laming et al. 2006; 
Hwang & Laming 2009). Cas A particularly lends itself to quantitative analysis by virtue of 
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its well-constrained age (Thorstensen et al. 2001; Fesen et al. 2006) and distance (Reed et al. 
1995). Recent one- dimensional hydrodynamical models for the interaction of the remnant 
with a circnmstcllar wind medium indicate that the reverse shock in Cas A has already 
interacted with a significant fraction of the ejecta (Chevalier & Oishi 2003). Laming & 
Hwang (2003) apply their models directly to Chandra X-ray spectra, and also infer that 
there is very little of the unshocked ejecta remaining. Thus a large fraction of the ejecta 
in Cas A are accessible to X-ray observations, the exceptions being the small amount of 
ejecta that have not yet been shocked, any ejecta that might have been subject to rapid 
cooling by thermal instability, and ejecta of very high density that have been shocked to 
lower temperatures and thus emit at longer wavelengths as optical knots (e.g. Fesen 2001). 

The organization of the paper is as follows: Section 2 begins by giving an overview of 
the analysis procedure, with explanations as to the spectral grid and background, and most 
importantly how the spectra were organized according to the fitting model required. We then 
proceed to discuss the results, from maps for the remnant as a whole, to a few observations 
about properties of the forward shock-dominated regions, and then the ejecta. Section 3 
discusses issues concerning the models and the refinements that have been made for this 
work compared to previous work. Section 4 then presents the results of the mass calculation 
and issues connected to it, such as the presence of unshocked ejecta and implications for the 
neutron star kick. 


2. Analysis Procedure 

We use the 1 Ms Advanced CCD Imaging Spectrometer (ACIS) observation obtained by 
the Chandra X-ray Observatory in February and May 2004 (one of the first Chandra VLP, or 
Very Long Project, observations). This observation is described by Hwang et al. (2004). Its 
most salient features are that it was obtained in nine observation segments (OBSIDS) and in 
GRADED mode, wherein the CCD detection events are characterized onboard the spacecraft 
before telemetry. It was necessary to use GRADED mode to reduce the telemetry load, due 
to the high source count rate (the observation accumulated some 280 million photon events 
in 980 ks). The main disadvantage of using GRADED mode is that detailed corrections 
of detector problems such as pulse pile up are not possible. The data were processed with 
CALDB Version 3.2.2. 

In this section, we outline the analysis procedure, which includes definition of the spec- 
tral grid, considerations for the subtraction of the background, the classification of forward 
shock and ejecta-dominated spectra, the selection of appropriate spectral models, and the 
results of the fitting. 



2.1. Spectral Grid 


To carry out our spectral survey of the ejecta in Cas A, we define a grid of 10" boxes 
across the entire supernova remnant (SNR), optionally subdividing by factors of two down 
to 2.5" depending on the number of counts. The regions are thus square or rectangular with 
sides of length 2.5, 5, or 10". We excluded the regions dominated by the central compact 
object, and a few other regions with low numbers of counts. This gave us 6202 spectra for 
analysis. Although this is a very large number, the data would have readily allowed at least a 
factor of two finer grid sizes down to nearly the Chandra angular resolution in some regions. 
This being the first attempt at an analysis of the complex thermal emission on such a fine 
angular scale, we tried to keep the number of grid regions relatively tractable. 

Each spectrum was then extracted individually from all nine OBSIDS comprising the Ms 
observation, and the corresponding individual photon redistribution matrices and effective 
area hies computed. The individual spectra were then added, and the individual response 
hies weighted and added, to obtain the hnal spectral and response hies for each grid region. 
The pulse-height spectra were further binned to yield a minimum of 25 counts per bin to 
apply Gaussian statistics. 

The filamentary and knotty characteristics of Cas A’s X-ray emission, its high surface 
brightness, and the excellence of the Chandra mirrors, result in pulse pile-up at the brightest 
ejecta knots. Pulse pile-up can produce spurious line features, for example, with energies 
that correspond roughly to the summed energy of the Si and/or S He a blends that are so 
prominent in Cas A. Pile-up of continuum photons will also skew the spectral shape to yield 
higher fitted temperatures. 

A detailed correction for pulse pile-up cannot be performed for these GRADED mode 
data, and we made no attempt to account for these effects for the thousands of spectra 
examined here. We can, however, estimate qualitatively the number of our spectra for which 
pile-up is a significant problem. For the vpshock fits described later (in section 2.3), we 
consider the distribution of the fitted temperature against the counts per pixel for each 
spectrum. The average fitted temperature shows a clear trend of increasing for spectra 
with more than about 6000 counts/pixel in 980 ks 3 . For the 25 regions with the highest 
count/pixel ratios (> 6000 counts/pixel), the average fitted temperature is 2.25 keV, which 
is significantly higher than the average of 1.7 keV for the remaining spectra. Inspection of 


3 The count/pixel cut-off we discuss here is of course specific for this particular observation. For reference, 
the brightest pixels contain on the order of 16000 counts, corresponding to a count rate of 0.05 counts per 3.24 
s frame. This corresponds roughly to pileup fractions of less than a few percent according to the Chandra 
Proposers’ Observatory Guide. 



6 


these spectra show the problems in fitting the line emission described above; these problems 
are not apparent for lower count-rate spectra. Since the pile-up occurs noticeably in only a 
limited number of regions (~25), we expect that its uncorrected effect on the total mass is 
not large. 


2.2. Background 

The off-source regions outside the SNR show a line-rich spectrum similar to that of 
the SNR itself. The background is thus dominated by the bright source spectrum due to 
scattering by intervening dust, the CCD readout, and to a lesser extent, the wings of the 
mirror point spread function. Dust scattering has been shown to be significant for Cas A 
based on ROSAT and Einstein observations by Predchl & Schmitt (1995) and Mauche & 
Gorenstein (1989). 

Since the source spectrum, and hence the background, varies both in form and brightness 
around the remnant, we extract 16 separate background spectra by subdividing a 0.4' thick, 
~ 3.5' radius circular shell surrounding the remnant so that each segment covers roughly 
20 degrees in azimuth. We then assign a background spectrum for each source spectrum 
based on the source region’s azimuthal position angle. Where the source extends beyond 
the background extraction radius, near the northeast “jet” and its southwest counterpart, 
we use the sum of the two nearest background regions to either side. 

We choose to subtract the background from the source spectrum rather than to fit a 
model to each background spectrum. A model for the background would have to include a 
thermal component for the scattered source spectrum, whose normalization would inevitably 
be difficult to constrain. The normalization for the particle background component is well- 
constrained by the spectrum at high energies, but the effect of scaling the scattered thermal 
background in the same way as the particle background is similar to simply subtracting the 
total background. 

Simple subtraction of the off-source background means that the background is generally 
underestimated, and by different amounts depending on the photon energy and the dominant 
scattering process. The precise determination of the background is thus a sophisticated 
problem that is beyond the scope of our study. The ejecta survey which is the focus of our 
efforts will be less sensitive to the precise background subtraction because most of the ejecta 
regions are bright, with strong line emission; the subtleties of the background subtraction 
are less important here than for inferences about the origin of weak lines. 
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2.3. Spectral Classification 

2.3.1. Spectral Map 

To provide a consistent spectral characterization of the entire remnant, we fit to every 
spectrum in our grid a single-component plane-parallel shock model with variable element 
abundances (XSPEC model trp shock: Borkowski et al. 2001) modified by interstellar absorp- 
tion (XSPEC model wabs ; Morrison & McCammon, 1983). The vp shock model incorporates 
ionization ages from zero up to a fitted maximum value. In previous work (Laming & Hwang 
2003; Hwang & Laming 2003), we had used simple nonequilibrium ionization (NEI) models 
characterized by a single temperature and ionization age, but there the extraction regions 
were generally as small as 3" and the exposure time only 50 ks. We find that for the much 
longer exposure time of these data, and the sometimes larger extraction regions, such models 
are no longer acceptable. Because most of the regions have ejecta-dominated spectra, we 
have not included any elements lighter than O in the model, taking O to be the primary 
source of the continuum in the manner of Vink et al. (1996), Laming & Hwang (2003) and 
Hwang & Laming (2003). The abundances of elements heavier than O are generally varied 
(with Ni linked to Fe, Ca to Ar, and sometimes Ar to S), while O is held fixed at the solar 
value, using the solar abundances of Anders & Grevesse (1989). The redshift is allowed to 
vary to either positive or negative values, given that the ejecta in Cas A have significant bulk 
motions (e.g., Markert et al. 1983; Willingale et al. 2002), and is generally driven by the Fe L 
emission and the strong Si blend. Given that there may be calibration problems near Si (see 
the discussion in DeLaney et al. 2010), we consider the fitting of the redshift more as an aid 
to achieving good fits rather than a reliable measure of the actual line-of-sight velocity. The 
fitted redshift does, however, reproduce the gross features of the bulk motion as measured 
by other means. A Gaussian smoothing scale for the spectrum was also fitted, given that 
the observed strong spectral lines are generally broad. 

While we do not present errors for the fitted parameters here, these were computed for 
key parameters. We have found that this is a necessary step to optimize finding the true 
minima: the \ 2 terrain is generally rather rugged for these fits. The redshift and ionization 
age in particular are prone to settle at secondary minima, as might be expected, considering 
that these parameters both have a strong effect on the energies of the strong line features. 

Here we note some of the general features of these results based on images shown in 
Figure 1, where fitted model parameters such as temperature, ionization age, and various 
element abundances are shown for each spectral region, and histograms and a scatter plot 
for the fitted temperatures and ionization ages shown in Figure 2. The column density 
increases systematically to the west by about a factor of 2 or more and reaches a maximum 



at the western extremity of the remnant. This is entirely consistent with previous work 
(Keohane et al. 1996; Willingale et al. 2002). There are also smaller regions in the center 
with significant localized column density enhancements. The temperature distribution has 
two peaks, seen most clearly in the histogram in Figure 2, with the lower peak associated 
with ejecta-dominated regions and the higher with FS-dominated regions. The ionization 
age, by contrast, is remarkably narrowly distributed overall, though coherent regions with 
high ionization ages can be identified on the map in Figure 1. In the western half of the 
remnant, these high ionization regions appear to be associated with the FS. The highest 
ionization ages are found in the east, however, and there they are clearly associated with 
Fe-enriched ejecta, as is apparent from comparing the ionization age and Fe abundance maps 
shown in Figure 1. 

We will focus our attention primarily on the abundance maps for Fe and Si, although 
such maps are also shown for all the elements for which we fitted the abundances. The 
maps for S and Ar are on the whole similar to that of Si. The abundance maps for Si and 
Fe are distinctly different, but both feature three main lobes of ejecta located north, east, 
and west. The Si is also extended to the northeast along the northeast ejecta “jet”, while 
Fe is particularly distinct in the eastern region compared to the Si, and shows its highest 
enhancement in the outermost parts of this region. In general, the abundance maps we 
obtain are strongly reminiscent of the corresponding line images shown by Hwang, Holt, & 
Petre (2000) and Hwang et al. (2004), and also resemble the abundance maps obtained by 
Willingale et al. (2002). Those authors have previously noted that the line emission and 
element abundance patterns for Si, S, Ar, and Ca in Cas A are similar to each other. 

The Ne and Mg maps show a distinct character in that they do not show any prominent 
morphological characteristics aside from a brightening at the western end of the remnant. 
One must be cautious to interpret those results, however, as the the line emissions of Ne and 
Mg fall in complicated parts of the spectrum, and the abundances of these elements may 
be correlated in complicated ways with parameters such as the ionization age, which shows 
similar distribution patterns in the west, and the column density, which is very high in that 
region. It is clear, however, that Ne and Mg both show a strikingly different morphology to 
Si, S, and Ar, or to Fe, and are much more similar to each other than to any of the other 
elements. 


2.3.2. Forward Shocked Regions 

Broadly speaking, the distinctions between ejecta- and FS-dominated regions in Cas 
A are readily apparent, with differences in temperature, ionization age, and element abun- 
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Fig. 1. — Partial set of fitted parameters for the single vpshock model fits. Top row: Maps of 
X 2 per degree of freedom (left), and of fitted column density Nh in 10 21 cm 2 (right). Middle 
row: Maps of temperature kT in keV and ionization age n e t in cm _3 s (sqrt scale). Bottom 
row: Maps of fitted Si and Fe element abundances, relative to the solar values of Anders & 
Grevesse (1989) by number, both on sqrt intensity scale, truncated at the high end. 
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Fig. 1. — Fitted abundances of the remaining elements, Ne, Mg, S, Ar. The intensities are 
square-root scale and truncated at the high end for S and Ar. 
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dances. To carry out a survey of the ejecta mass, however, we must either model both the 
ejecta and CSM components for each spectrum or else identify the specific regions where 
the reverse-shocked ejecta make the dominant contribution to the emission. Given the scope 
of the spectral analysis, we have adopted the latter approach. Multi-component fits can be 
difficult to constrain reliably, particularly if one of the components is relatively weak, and 
thus would require more individual attention than is feasible for a sample of thousands. Con- 
sequently, our next aim is to identify and eliminate regions whose spectra can be completely 
associated with the forward shock. 

We evaluate the presence of thermal emission associated with the forward shock by 
fitting a second set of plane-parallel shock models to every spectrum, but this time with 
element abundances appropriate for the CSM. The optically emitting quasi-stationary floc- 
culi (QSFs) in Cas A are understood to be circumstellar mass loss from the progenitor, but 
abundance measurements for QSFs are limited to a small number of knots. These show 
an order of magnitude enhancement of N and sometimes also of He (Chevalier & Kirshner 
1978). Theoretical calculations for the presupernova composition are also given by Arnett 
(1996), where the models allow the elements H, He, and N, all apparently present in Cas 
A, to exist simultaneously at a narrow temperature range near log T (10 9 K) = 1.5 (their 
Figure 7.6). At that temperature the abundance of He is 3 times the solar value of H by 
number, and that of N about 15 times the solar value. As these abundances for He and N 
are broadly consistent with the observational measurements, we proceed to adopt them for 
our fits, along with solar values for the remaining elements, as representative CSM element 
abundances. 

About 1209 regions gave reasonably good fitting results (y 2 < 1.2) with the vpshock 
model and these element abundances, and are thus assigned to the forward shock. They are 
distributed mainly in the remnant’s outer rim and southwest interior, as would be expected 
based on the 4-6 keV X-ray continuum image that highlights the forward-shocked regions 
(Gotthclf et al. 2001). Their average temperature is 2.2 keV, and their ionization ages 
are rather narrowly distributed with an average value of 2 x 10 11 cm _3 s. These values 
can be assessed in the context of the models of Laming & Hwang (2003), which give the 
current density of the CSM at the forward shock at about 1.5-2 cm -3 . Considering the r~ 2 
dependence of the circumstellar density, the forward shock will have encountered much denser 
material in the past and the present-day ionization state of the forward shocked material 
is expected to be relatively advanced. The models give values of the ionization age in the 
10 11 cm -3 s range, approaching 10 12 cm _3 s; they also indicate that gas is also rather hot, 
with temperatures from 2.5-4 keV. The average values of the temperature and ionization age 
that we find in our regions are close to the range predicted by the models, though the fitted 
spectra don’t show as broad a range in ionization age as is expected. 
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To the forward-shock regions identified solely by the thermal emission model above, we 
must also add those that have a strong nonthermal contribution. To identify these, we devise 
a rough diagnostic for the smoothness of the X-ray spectrum. We bin each background- 
subtracted spectrum at each significant line feature and continuum interval (some of these 
cover only a narrow energy range), compute the ratio of counts for each major line feature 
relative to counts in an adjacent continuum bin, and take the sum of these ratios for all 
the line features. The distribution of this quantity has two overlapping peaks; we take the 
spectra associated with the lower peak (corresponding to weak lines in the spectrum) and 
perform further fits with a composite plane-parallel pshock plus power-law model. For a 
cutoff in y 2 < 1.2 for these fits, we associate 206 additional regions with the forward shock, 
giving a total of 1415 regions with spectra that are consistent with emission associated with 
the forward shock alone. 

Figure 3 summarizes the distribution and spectral characteristics of the forward shock 
regions. Their locations echo the 4-6 keV continuum maps shown by Gotthclf et al. (2001) 
and Hwang et al. (2004). As already noted, their temperatures and most especially their 
ionization ages are rather narrowly distributed, more so than for the sample as a whole. 
From the Figure, the peaks in the distributions correspond roughly to kT = 2.2 keV and log 
(n e t) = 11.25, or n e t = 1.8ell cm _3 s. These regions will not be considered further here as 
we focus our study on the ejecta. We will undertake a detailed consideration of the thermal 
emission associated with the forward shock in subsequent work. 


2.4. Spectral Survey of the Ejecta and the Presence of Pure Fe 

We associate with ejecta the remaining, more than 4000, regions that are inconsistent 
with forward shocked material. For these regions, the basic spectral model is the simple 
one-component plane-parallel shock that has already been presented and discussed. As 
noted above, we take the view that, for the ejecta sample as a whole, it is justifiable to 
assume sufficient ejecta dominance to neglect the forward shock component. Representative 
examples of all the various types of spectra seen in Cas A are shown in Figure 5. 

For some ejecta regions, the one-component spectral model is clearly inadequate to 
describe the ejecta emission in that it fails to account for the Fe K blend. This has already 
been noted by Hwang & Laming (2009). While there may be significant deficiencies in the 
atomic data that are used for the thermal emission models, particularly in the nonequilibrium 
ionization case, these cannot explain so large an effect as we see. In some cases, a strong 
Fe K blend is completely unaccounted for by the model, which otherwise characterizes the 
spectrum well. 
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Fig. 2. — The distribution of temperature kT and ionization age n e t amongst the 6202 
spectral regions fitted with pshock models (see Figure 1) : (left) two-dimensional distribution; 
(middle) histogram for temperature kT ; (right) histogram for ionization age n e t 



Fig. 3. — (left) Contours showing the 1415 regions that have been associated with the forward 
shock, and hence excluded for the ejecta mass calculation, based on the merged QSF and 
QSF+PL fits. For these forward shock regions is shown the distribution of fitted (middle) 
temperature kT (the scale is truncated above 10 keV), and (right) log 10 ionization age n e t. 
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Given the strong chemical inhomogeneities seen and expected in core-collapse supernova 
remnants, it is plausible that Fe-rich ejecta could be superposed along a given line of sight 
with ejecta of more normal composition, and that each would have distinct plasma conditions 
as well as abundances. Moreover, nucleosynthesis models predict the formation of nearly 
pure Fe ejecta by a-rich freeze-out during complete Si burning; this is in addition to ejecta 
composed of a range of elements from Si and Fe formed by incomplete Si-burning. 

Hwang & Laming (2003) and Hwang et al. (2004) have reported the presence of highly 
enriched Fe ejecta in the southeastern region of Cas A, with one region particularly appearing 
to be a “pure Fe” cloud. We take another look at this region using the 1 Ms observation. 
Figure 4 shows a portion of the southeastern region of Cas A and the roughly 2.5" extraction 
region that we use for this cloud, which is smaller than that used by Hwang et al. 2003. The 
background is taken locally from the vicinity of the source region, and we compared results 
using three different background spectra. The background was modeled in this case, rather 
than subtracted, and the model included in the fit to the unbinned, unsubtracted source 
spectrum. The results for all three background subtractions gave fitted Fe/Si abundance 
ratios that are enhanced over the solar values by factors of 10-20 by number (or factors of 
20-40 by mass). The right panel of Figure 4 shows the fit giving the highest Fe/Si ratio. 
Note the absence of recognizable emission from S or Ar, and the weakness of the Si emission; 
the emission lines at energies above 7 keV are of H-likc Fe. These results demonstrate that 
it is possible to identify emission in Cas A that is plausibly associated with pure Fe. 

For a subset of 700 regions with the strongest Fe line emission and the poorest fits in 
the Fe K band, we tested two options for the added component representing either shocked 
circumstellar medium or pure Fe and Ni ejecta. For the former, we take a plane-parallel 
shock with CSM abundances and the average temperature (2.2 keV) and ionization age 
(3 x 10 11 cm _3 s) obtained for the forward shock regions, all as described above. For the 
latter, we found that a simple nonequilibrium ionization component, with a single value 
of the ionization age was sufficient, with the temperature fixed at 1.95 keV and the single 
value of the ionization age fixed at 8 x 10 11 cm _3 s. These values are representative of the 
most enriched Fe ejecta spectra, such as those studied by Hwang & Laming (2003); the fits, 
however, are not strongly sensitive to the exact values used so long so long as they are of this 
approximate magnitude. We fix the parameters of the second component whenever possible 
because the limited spectral resolution makes it very difficult, even at these exceptionally 
high levels of signal-to- noise, to constrain both model components independently. As before, 
errors were calculated for key parameters to help ensure finding the true minimum in y 2 . 
Only one of these 700 ejecta+CSM fits gave y 2 < 1 in the 6-7 keV Fe K region, leading 
us to conclude that the presence of a forward shock component can not explain the strong 
Fe K emission. On the other hand, a majority of the vpshock+NEl fits did give acceptable 
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Fig. 4. — (Top) Southeast region of Cas A showing the extraction region used here for the 
pure Fe cloud of Hwang & Laming (2003). (Bottom) Fitted spectrum, in which the local 
background has been modelled rather than subtracted and the unbinned spectrum fitted 
using C-statistics. The spectrum is shown binned for clarity, with the black trace giving the 
total source plus background spectrum, and the red trace giving only the source spectrum. 
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X 2 values in the 6-7 band, although the improvement was sometimes dne to a better match 
to the continuum rather than the line emission. An example is shown in final two panel 
rows of Figure 5, where the same spectrum is shown with the original single-component 
vpshock fit, a two component vpshock fit for normal ejecta and shocked CSM, and finally a 
two ejecta-component fit where the added ejecta component is taken to be pure Fe (and Ni). 

On a broader scale, we identify 3032 spectra where a composite ejecta model may be 
needed by selecting regions where the fitted Fe abundance for the basic single-component 
model was above 0.3 solar, and focus only on two component models with the usual vpshock 
ejecta component and a second (NEI) ejecta component corresponding to Fe. We evaluate 
using an f-test those regions where the improvement in fit from adding the second component 
was above a certain probability threshold. The calculation of the mass of various elements 
in the ejecta is discussed fully in Section 4. Here, we simply show in Figure 6, examples 
of the Fe mass distribution in the second, “pure” Fe component for probability cuts of 0.1, 
0.01, and 0.001 for adding this component. The spatial distribution retains its main features 
regardless of the exact probability cut imposed. 


3. SNR Models 

In this section we provide an overview of our models for the evolution of Cas A. These 
are treated in a manner similar to that in Laming & Hwang (2003) and Hwang & Laming 
(2003), for expansion of the remnant into a circnmstellar wind (density oc r~ 2 ), but with 
modifications that we discuss here. 


3.1. Dynamics 

We modify the models to accommodate SNR expansion into a stellar wind “bubble” 
as presented and discussed in Hwang & Laming (2009). The chief motivation involves the 
interpretation of infra-red light echoes associated with Cas A as reprocessed shock breakout 
radiation (Dwek & Arendt 2008). This radiation is concentrated in the UV-ELTV spectral 
region, and would be strongly absorbed in photoionizing any intervening neutral gas, which 
presumably would be a relic of the red supergiant presnpernova stellar wind. Most of this 
opacity would lie close to the progenitor, and so a small “bubble” in the circnmstellar medium 
would allow the breakout radiation to escape and illuminate the surrounding dust. 
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Fig. 5. — Sample spectra exemplifying typical spectral types seen in Cas A (first two rows) 
Spectra associated with the forward shock, showing various mixtures of thermal emission 
associated with shocked circumstellar medium and nonthermal emission, (third row) Typical 
“Si-dominated” and Fe-dominated ejecta spectra, (last two rows) For the same spectrum, 
a comparison of single-component vpshock ejecta models, ejecta plus shocked CSM models, 
and a model with two ejecta components corresponding to “normal” ejecta and “pure Fe” 
ejecta. 
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Fig. 6. — The distribution of the added “pure” Fe ejecta component for f-test probability 
cutoffs of 0.10, 0.01, and 0.001 (left to right). In the left panel, we show fiducial circles 
showing the radius and position of the forward shock (2.6 pc at 3.4 kpc distance, in yellow) 
and the contact discontinuity (1.9 pc, in white), and the velocity vectors for the NS (white), 
the “pure” Fe ejecta (light blue), and total ejecta (yellow) from Table 2 (see subsection 4.2). 
The circles and velocity vectors are centered at the explosion center of Thorstensen et al. 
( 2001 ). 
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Another update corrects the motion of the reverse shock 4 . In previous work we have used 
a model of the ejecta where a uniform density core is surrounded by a power law envelope. 
The evolution while the reverse shock is in the outer envelope is given by Chevalier (1982); 
this solution is then coupled to the known asymptotic behavior of the Sedov-Taylor solution. 
In this phase, we have previously held the reverse shock speed constant at its value at the 
envelope-core transition, although it is known that some small acceleration should occur 
during the core propagation. We now implement a reverse shock acceleration in this phase 
of evolution following the results given for a suite of models by Patnaude & Fesen (2009). 
By matching to their results, we write for the reverse shock speed and radius: 


v r 

R r 


— 3-s i’b I n — 6 j. 
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where v r = R r /t — dR r /dt. Here s is the power law index for the CSM density pro- 
file, n that for the ejecta envelope, Rb and Vb the forward shock radius and speed, re- 
spectively, l ED the “lead factor” (the ratio of forward to reverse shock radii, R b /R r ), and 
Xq = (40 .74:M e j/ pRl) 1 ^ 3 -^, with p being the density of the CSM at the forward shock in H 
atoms (or equivalent mass) per cm 3 , t core is the time in years when the reverse shock hits 
the ejecta core, and M e j is the ejecta mass in solar masses. 


In Table 1 we give various models in the range of plausible ejecta masses 2-4 M 0 , all 
designed to match as far as possible the known dynamics of Cas A. We take a distance of 
3.4 kpc (Reed et al. 1995) and an explosion date of 1671.3 ± 0.9 as determined from the 
proper motions of 17 high velocity outer ejecta knots (Thorstensen et al. 2001). These knots 
are identified on plates taken as long ago as 1951 and are assumed to be undecelerated, 
providing the longest time baseline for such studies. More recently, Fesen et al. (2006) have 
determined a best explosion date of 1681 ± 19 from a sample of 126 relatively undecelerated 
knots located along the NW limb, observed over a time baseline of 9 months from March to 
December 2004. With 1671 as the earliest possible explosion date, the remnant age may be 
up to 333 years in 2004. 


Further constraints on the model are provided by the observed shock velocities and radii. 
A number of authors have made measurements of the expansion rate of Cas A, both for its 
bright ejecta ring (Vink et al. 1998; Koralesky et al. 1998; DeLaney et al. 2004) and the 
forward shock (DeLaney & Rudnick 2003). The most recent summary is given by Patnaude 


4 Some other minor modifications and corrections to Appendix A in Laming & Hwang (2003) are collected 
in the appendix to this paper 
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& Fesen (2009), who incorporate measurements over the longest time baseline (2000-2007). 
They find that the expansion is slightly slower in the N and NW filaments, consistent with 
the known presence of denser CSM on the western limb of Cas A. We take the average 
expansion 0.31” yr -1 from other parts of the remnant, which (at a distance of 3.4 kpc) gives 
a forward shock velocity of 5050 km s -1 . Including all measurements gives a slightly lower 
average expansion of 0.30” yr -1 , and an inferred forward shock velocity of 4850 km s -1 . 

The forward shock radius is measured from the wispy structures seen in 4-6 keV contin- 
uum emission. Gotthclf et al. (2001) use data taken in 2000 and determine an average radius 
for the northwest sector of 153 ± 12”, which translates to 2.52 ± 0.2 pc, but the true shock 
radius should probably be determined from the outermost of these of these structures rather 
than the average. Helder & Vink (2008) perform a deconvolution for various azimuthal 
angle ranges and give a slightly larger forward shock radius of 160” (2.637 pc). There is also 
a small variation in the radius around the limb. We take a value of 2.6 pc for the forward 
shock radius in 2004 when the Chandra VLP observation was taken. 

The reverse shock radius is harder to quantify. Gotthelf et al. (2001) and Helder & Vink 
(2008) give a range of radii 1.52 - 1.73 pc, depending on location in the remnant. Morse et 
al. (2004) measure the reverse shock velocity from proper motions derived from two Hubble 
Space Telescope WFPC2 images separated by two years. At a reverse shock radius of 106” 
- 122” (1.75 - 2.01 pc) the reverse shock velocity with respect to the expanding unshocked 
ejecta is of order ~ 2000 km s -1 . 

The close match of Cas A’s forward shock expansion parameter of 0.66 to the Sedov- 
Taylor value for expansion into an inverse square density profile has focussed interest in low 
ejecta mass models that allow faster evolution to the Sedov-Taylor limit. Even so, models 
generally predict too high a shock velocity and/or too small a radius for the forward shock 
compared to observations. To increase the forward shock radius, it has been suggested 
(Hwang & Laming 2009, and see above) that Cas A exploded into a small bubble in the 
circumstellar medium that developed with a transition from a slow dense red supergiant 
wind to a fast tenuous wind shortly before explosion. Another possibility is that cosmic 
ray energy losses at the forward shock provide extra deceleration to suffciently reduce the 
expansion parameter (Patnaude & Fesen 2009). 

The plasma density ahead of the forward shock is derived by Willingale et al. (2003) 
from fits to XMM data. The baryon mass and filling factor quoted in their Table 3 give 
a preshock density of hydrogen atoms or equivalent mass of 1.47 cm -3 , assuming a shock 
compression of a factor of 4. Their quoted value for the mass of shocked CSM, 8.31 M 0 , 
combined with the dimensions outlined above gives a larger preshock density of 1.99 cm -3 . 
These values give the time invariant quantity pr‘j = 10 — 13.5 cm _3 pc 2 where the density p 
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at the forward shock is in H atoms (or equivalent mass) cm 3 and the forward shock radius 
Tf is in pc. 

Similar results are given by the fits of thermal bremsstrahlung spectra to Suzaku data 
(Maeda et al. 2009), interpreted as emission from shocked circumstellar medium. The 
preshock ion number density derived from their thermal bremsstrahlung and power law 
fit is 2.1/ ( Z 2 ) cm -3 , where ( Z 2 ) = 'Yh i n i Z^ / ’Y^ i n i Z i is an average charge per ion coming 
from the thermal bremsstrahlung emissivity. For gas composed of equal number densities of 
H and He, (Z 2 ) = 1.7 and the preshock density is 1.24 cm -3 in H atoms of equivalent mass. 
Taking fit results using thermal bremsstrahlung and SRCut or a cut-off power law gives a 
higher density of 1.5 cm -3 and prj = 8.4 — 10 cm _3 pc 2 . 


3.2. Cosmic Ray Energy Losses 


The expansion parameter for a radiative blast wave can be different from the Sedov- 
Taylor limit in the case that radiative losses from the shock are significant, or if the shock 
dissipates energy via cosmic ray acceleration. The two cases may be treated in a similar 
manner. Following Liang & Keilty (2000) and generalizing to the case of a stellar wind 
preshock density profile p oc r~ s , the radius of a radiative blast varies with time as 


r oc 


( 3 ) 


where 7 is the adiabatic index of the gas, a 


7 + \/( 2 - if + 4 (71 - 1) | / 4, 


and 


71 = | a /1 + e (7 2 — 1) — 1 1 x 2/ (7 — 1) e — 1 is the adiabatic index modified by the loss 
of a fraction e of the shock energy to radiative (or cosmic ray) losses. If e = 0 so that 
7 = 71, then a = 1/2 for all 7 and r oc t 2 ^ 5 ~ s \ Note that energy going into cosmic rays 
that remain trapped with their energy in the shock does not change the dynamics, except 
through modifying 7. If all the shock energy is radiated (e = 1), then a = 1 — 7/2, and 
the expansion parameter changes from 0.4 to 0.28 (for s = 0) and from 0.67 to 0.55 for 
s = 2, taking 7 = 5/3. Consequently for more modest cosmic ray energy losses, rather small 
deviations from the unmodified expansion parameter are expected. Cosmic ray energy losses 
of 10% only change the s = 2 expansion parameter from 0.67 to 0.65, decreasing to 0.63 
for 30% losses. Patnaude & Fesen (2009) find variations in the expansion parameter of this 
order. 



24 


3.3. The Ejecta Filling Factor 

A correction is needed for unresolved density structure in the ejecta. Models of super- 
nova remnants expanding into a stellar wind (p oc r~ 2 ) density profile generally show a strong 
density “spike” in the radial direction at the contact discontinuity, and if the reverse shock 
is still propagating through the outer power law portion of the ejecta density profile, the 
density at the contact discontinuity is formally infinite (Chevalier 1982). Later in its evolu- 
tion, the density spike is still present (Chevalier & Oishi 2003), with a width corresponding 
to about 2 arcsec. This is comparable to or smaller than the size of many of our fit regions, 
so we expect that ejecta structures will not always be resolved, and a filling factor correction 
is therefore needed when calculating the plasma electron density (see subsection 4.1 below). 
We take the volume of each spectral region o be A 3//2 where A is the area imaged in arcsec 2 . 
This is modified to A 3//2 /, where / is the “filling factor” of the ejecta density profile. 

We estimate a quantitative filling factor / to be applied to our fit results in terms of the 
density scale length L p as / = L p /\/A. 5 In deriving L p , we follow Hamilton & Sarazin (1984), 
who give a similarity solution for ejecta structure close to the contact discontinuity. Treating 
the expansion of uniform density ejecta, they use a coordinate system with z = 0 at r = v exp t , 
the outer limit of the ejecta at time t if they were freely expanding without interacting with 
the reverse shock. The location of the reverse shock is z s , giving dimensionless Eulerian and 
Lagrangian coordinates ( = zjz s and Co = z o / z s respectively, where z 0 is the initial value of 
z in a Lagrangian plasma element at the time that the plasma was shocked, i.e. z s at the 
time of shock passage. With these definitions, close to the contact discontinuity they find 
(their equation 13) 



where v = 23/15 for a stellar wind external density profile, ( c is the value of C at the contact 
discontinuity, and V s = 1/12 is the dimensionless reverse shock velocity. The pressure is p, 
and p s is the pressure at the reverse shock. Neglecting the dependence of p on r in this 
region, we write 

dp _ dp d( 0 gC _ 6 p f p _\ 3/5 _ 6 p p 

dr d( 0 d( dr 5 z a V a \p s ) 0 5z s u (C - Cc) L p 


5 We are considering emission from a “sheet” of plasma, so only one dimension is compressed, rather than 
from a clumpy medium, where all three dimensions would require correction. 
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so that the density scale length L p = 23 — ( c ) z s /18. We calculate the average of this 

quantity (L p ) = JJ C pL f) r 2 dr / JJ C pr 2 dr using 


I pr 2 dr = jj; ( Ps(q 6/5 r 2 dr = p s M 
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The integrals are evaluated using the substitution u — r c — r, with the final result 

(1 - r s /r c ) 28/23 - f (1 - r s /r c ) 51/23 + f (1 - r s /r c ) 74/23 
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Including the variation of p, L p in equation 5 is modified to 


f 5z//3 — dlnp/dlnC'O 1 
| 2 — <91np/<91nCo J 


( 8 ) 


From Table 1 of Hamilton & Sarazin (1984), <91n/)/<91n((o varies from -1 at the reverse shock 
to about -0.6 at the contact discontinuity, which reduces the value for L p derived assuming 
<9 lnp/<9 In <(o = 0 by 5% to 7%. For completeness, if the external density oc r -3 , then 
v = 6/5 and L p = z s (( — ( c ) = r c — r, independent of <91np/<91n</ 0 . The filling factor is then 
calculated as / = L p /\/~A, with L p and A here expressed in arcsec and arcsec 2 respectively. 

As we will show below, our data suggest that many secondary shocks reflect from 
the blast wave and reverse shock as they encounter density inhomogeneities in Cas A. 
We therefore use equation 7 with the tacit assumption of constant pressure near the con- 
tact discontinuity. Truclove & McKee (1999) point out that this approximation is also 
made by Hamilton & Sarazin (1984) in their energy equation. With r s /r c = 0.9 we find 
L p = 0.021r c = 1.26 x 10 17 cm, which at a distance of 3.4 kpc subtends an angle of 2.48 
arcsec. In this last step we are extrapolating slightly out of the formal range of validity of 
the Hamilton & Sarazin (1984) model. There, z c / z s = (v exp t — r c ) / [y exp — r s t) = 0.953, 
which requires an expansion velocity of the outermost ejecta greater than 20,000 km s^ 1 
for values of r c = 1.95 pc and r s = 1.7 pc appropriate for Cas A. It is also worth bearing 
in mind other caveats to the application of this model. Hamilton & Sarazin (1984) assume 
uniform density ejecta, while we adopt the prescription of Truelove & McKee (1999) of a 
core-envelope density profile. We have also introduced the circumstellar “bubble” into which 
Cas A is assumed to have exploded. 
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4. Discussion 

4.1. The Elemental Composition of the Cas A Ejecta 

In previous works (Laming & Hwang 2003; Hwang & Laming 2003) we have discussed 
the use of models of SNR evolution to interpret the fitted ionization age of ejecta spectra 
in terms of the time elapsed since reverse shock passage, hence allowing the inference of 
a Lagrangian mass coordinate. This is successful in a limited number of carefully selected 
regions of Cas A, but is evidently not the case in general. Figure 2 shows the ionization age 
value strongly peaked near log (n e t) — 11.3, which implies that all the ejecta are piled up 
at a mass coordinate of 0.2 (where 0 marks the center of the ejecta and 1.0 the outermost 
extent), instead of being uniformly distributed with mass coordinate. We suggest that 
this must arise from the interaction of the reverse shocked ejecta with secondary shocks 
propagating within the SNR shell. These arise as the forward and reverse shocks encounter 
density inhomogeneities, sending a transmitted shock through the density structure together 
with a reflected shock back into the previously shocked gas. Such a scenario has previously 
been considered (Laming 2001a, b) as a possibility for electron acceleration with a view to 
explaining the hard X-ray emission of Cas A. More recently, Inoue et al. (2010) consider a 
similar model for ion acceleration at shocks. 

Taking a typical electron density of 200 cm -3 (Lazendic et al. 2006, and see below), an 
ionization age of 2 x 10 11 cm _3 s implies a shock interaction around 30 years prior to the 
observation. Given this “obscuration” of the ionization age determined by the reverse shock, 
we estimate the density in each fit region from the emission measure assuming 

EM = nin e ZfV = n e V r kZ 2 = n\V (Z 2 ) (9) 

i i 

where (Z 2 ) = )>T n^Z 2 / )>T 71,2) = JTnjZ 2 /n e , and the sum i is over all ions in the plasma. 
The plasma volume is V . The plasma electron density then follows with an appropriate 
assumption about V. The simplest assumption that V = A 3 / 2 (where A is the area on the 
sky of the fit region chosen), and we modify this to V = A 3 / 2 x / where / is the filling factor 
discussed previously. 

The ejecta mass in each imaged region is estimated as EM/n e / (Z 2 ) x ni r mi/n e , from 
which the masses of individual elements are obtained as fractions of the total ejecta mass 
from the fitted element abundances. In Table 2 we give the masses of O, Ne, Mg, Si, S, Ar, 
and Fe resulting from our fits, for assumed filling factors corresponding to a density spike 
thickness of 2.5” and 5”; the latter arises as one observes the SNR ejecta shell in both front 
and back along the same line of sight. In this latter case, the complete X-ray emitting ejecta 
mass inferred is 2.84 M &1 and the total ejecta mass including unshocked ejecta is 3.14 M & 
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with reference to Table 1; this gives a value of prj from Table 1 at the upper end of values 
inferred elsewhere as discussed in section 3.1. Extrapolation to other cases is simple, since 
element masses are proportional to y/J. Table 2 also gives plane of the sky velocities for each 
element (calculated from the position of the center of mass of each element, and assuming 
homologous expansion), and for the SNR as a whole, and compares abundance ratios relative 
to O with solar system values. 

We also comment here that Cas A is thought to have retained some ff at the time of 
explosion (Fesen 2001; Fesen & Becker 1991; Chevalier & Oishi 2003), and this may contribute 
to the thermal bremsstrahlung emission here attributed to O. In this case, our ejecta mass 
estimate would be an underestimate, and if the H were not uniformly distributed, the recoil 
velocity we attribute to the O ejecta could be in error. It is difficult to assess an uncertainty 
here, but we note that the amount of H involved is likely rather small and unlikely to alter 
our conclusion of an approximately 3 M 0 ejecta mass. 

In Table 3 we give the masses of Fe associated with incomplete Si-burning, and the 
“pure” component, the latter arising either from complete Si-burning or a-rich freeze out, 
for a variety of assumptions about the statistical significance of requiring the second fit 
component. We calculate the F-distribution following Press et al. (1992) and evaluate 
the probability that the second component is not justified by the improvement in y 2 , i.e. 
the probability that the residuals from the fit with and without the second component are 
statistically from the same distribution. Note that in Table 2 we have taken those Fe masses 
corresponding to a probability of 0.01. 


4.2. The Mass and Distribution of Fe 

Our analysis of more than 4300 ejecta regions gives a mass of Fe in the Cas A ejecta 
of ~ 0.08 — 0.14 Mq, depending on the assumed filling factor and threshold probability for 
acceptance of the “pure” Fe component. This mass is comparable to that expected (Eriksen 
et al. 2009), but gives only the Fe detected in the reverse shocked X-ray emitting ejecta. 
According to our models, an extra 0.18 - 0.3 M 0 of ejecta may be unshocked, interior to the 
reverse shock, and visible primarily in the infra-red. 

Infrared observations with Spitzer do show emission at the center of the remnant that is 
believed to be from unshocked ejecta, but of [Si II]. As discussed by DeLaney et al (2010), it 
is coincident with free-free absorption seen in the radio by Kassim et al. (1995), and appears 
to correspond to cool (< 1000 K), low density gas that has been photoionized. DeLaney et 
al. conclude, however, that it is unlikely that a significant fraction of these unshocked ejecta 



could be Fe. Relatively little infra-red Fe emission is observed at all in the remnant, and 
the main [Fe II] line at 26 // is blended with [O IV]. Isensee et ah (2010) demonstrate that 
[O IV] is more plausibly the dominant component of this blend, though they do claim a 2 o 
detection of [Fe II] from the center when all spatial bins are summed. Unblended [Fe II] 
emission at 17.94 /mi (Ennis et al. 2006) and 1.64 /mi (Rho et al. 2003) is also detected, but 
appears to be associated primarily with the bright ejecta ring, and not the interior. [Fe II] 
lines are also absent or very weak in optical and near infra-red observations of Cas A ejecta 
knots. Hurford & Fesen (1996) only report marginal detections of [Fe II] 8617 A in spectra 
of fast moving knots (FMKs), and Gerardy & Fesen (2001) report similarly on the [Fe II] 
lines between 1 and 2 /mi in FMKs, although these transitions are easily detected in the 
spectra of quasi-stationary flocculi (QSFs). Eriksen et al. (2009) do detect [Fe II] 1.257 /mi 
and 1.644 /mi lines from 15 FMKs and 4 QSFs, from which they estimate the reddening. 

Further consideration of the dust in Cas A sheds light on the possible presence of 
unshocked Fe. Large masses of cold dust had been reported in Cas A from submillimeter 
observations (Dunne et al. 2003), but it has been pointed out that some of the detected 
emission originates from foreground molecular clouds (Krause et al. 2004; Wilson & Batrla 
2005). This emission can also be explained with < 10^ 3 M 0 of conducting dust needles 
(i.e., Fe) that are formed in the ejecta (Dwek 2004). Rho et al. (2008) favor FeO dust to 
explain the Spitzer infrared spectra and conclude that Fe dust may be present at masses 
up to 1O~ 2 M 0 . More recently, Nozawa et al. (2010) suggest the presence of significant 
quantities of cool dust, a conclusion that is supported by recent infrared observations with 
Herschel (Barlow et al. 2010), AKARI and BLAST (Sibthorpe et al. 2010). The cool dust 
has a temperature of about 35 K based on the infrared flux densities, and is presumably 
unshocked ejecta as it is confined to the central regions of the remnant. The emission is 
consistent with a silicate dust composition, and the inferred dust mass is about 0.06-0.075 
M 0 . Nozawa et al. (2010)’s calculations are for dust formation in a Type lib event (with 
an eye toward Cas A) and indicate that little dust would be associated with the innermost 
Fe-Ni layer due to the extended radioactive heating: the gas density drops too low before 
the temperatures are low enough for Fe or Ni grains to condense. Lip to about 10 ” 3 M 0 of 
Fe could be locked up in FeS grains, but highly Fe-rich ejecta would not necessarily have a 
dust signature, as appears to be borne out by the infrared observations. Cherchneff & Dwek 
(2010) take a chemical kinetic approach and find rather more FeS grain formation in a 20 
M q model SN with unmixed ejecta (the most appropriate of the cases they consider for Cas 
A), with about 0.021 M 0 of Fe in FeS. A mass this high is unlikely for Cas A, however, given 
that its highly stripped progenitor underwent a Type lib event. Dust will form relatively 
more efficiently in a normal Type II event because the overlying stellar envelope restricts 
ejecta expansion and preserves higher densities favorable for dust condensation. 
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We are left with the surprising conclusion that almost all of the Fe ejected by the 
supernova is now well outside the reverse shock and visible in X-rays, with very little left in 
the center of the remnant. This is true even for the “pure” Fe component, which we take to 
be the ashes of ct-rich freezeout. An initial mass of 44 Ti of 1.6 x 1O _4 M 0 has been inferred for 
Cas A from observations with the Compton Gamma Ray Observatory/COMPTEL (Iyudin 
et al. 1994), BeppoSAX and the INTEGRAL IBIS/ISGRI instrument (Vink et al. 2001; 
Renaud et al. 2006). Based on spherically symmetric simulations, designed to match the 
Cas A progenitor (Young et al. 2006; Eriksen et al. 2009), Magkotsios et al. (2010) predict 
~ 1O“ 4 M 0 of 44 Ti, and 0.25 M 0 of 56 Ni. These are both within a factor of 1.6 of measured 
values (taking the 5” filling factor in Table 2), but are discrepant in opposite directions so 
that the mass ratio 56 Ni/ 44 Ti~ 2500 compared with the observed ratio of ~ 10 3 . A similar 
outcome of the earlier calculations of The et al. (1998) and others (reviewed in The et al. 
2006) led Nagataki et al. (1998) to suggest that this ratio may be reduced in an asymmetrical 
explosion. Magkotsios et al. (2010) also enable an estimate of the Fe formed in ct-rich freeze 
out compared with that formed in incomplete Si- burning. Their Figure 17 suggests that 
approximately half the Fe should be associated with a mass ratio 56 Fe/ 28 Si> 10 and half 
with 56 Fe/ 28 Si< 10. In our Table 3, taking probabilities between 0.1 and 0.01, we would 
infer that approximately 1/3 of the Fe should have formed in a-rich freeze out and thus be 
in regions with very small 28 Si masses. 

In Cas A, the ionization age of > 10 11 ' 3 cm -3 s places the “pure” Fe at least 0.75 M 0 
out from the center in 3 M 0 of ejecta. According to Magkotsios et al. (2010), this would 
place it well outside the ejecta region where it was formed, into the outer layers where the 
composition would originally have been dominated by O, Si, C, and Ne. The fact that 
essentially all of the Fe is found at such locations, and not in the center of the remnant 
inside the reverse shock suggests the operation of a strong instability. 

The 44 Ti should also be expected to be well outside the reverse shock, spatially coincident 
with the pure Fe, which should be confirmed by observations with NUSTAR (Harrison et al. 
2010). This is quite different to the case of SN 1987A, where the ejecta are currently heated 
by the decay chain of 44 Ti, which appears to be located in the center of the remnant, with 
ballistic expansion velocity of order 3000 km s _1 (Kjaer et al. 2010). Gl.9+0.3 (Borkowski 
et al. 2010) is more like Cas A in this respect, with most of its 44 Ti in an X-ray bright region 
beyond the reverse shock, and a smaller mass of 44 Ti in the unshocked interior. Considering 
the young age and fast expansion of Gl.9+0.3, this may mean that 44 Ti has been mixed out 
even further into the ejecta than is the case in Cas A, although the SN origin of Gl.9+0.3 
is unknown. 

We are unable to derive reliable line-of-sight velocities from the Chandra data, but some 
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insight into the motion of the “pure” Fe component can be gained from the motion of 44 Ti 
observed with INTEGRAL (Martin & Vink 2008; Martin et al. 2009). The SPI instrument 
measures the width and the shift of the 1157 keV decay line of 44 Ca, which is the final step 
in the chain of decays from 44 Ti. It is found that the 44 Ti is receding from the observer with 
a velocity of about 500 km s” 1 , and has an intrinsic width, or velocity dispersion, of 430 
km s” 1 . Since the expansion velocity of the Fe is close to 4000 km s” 1 , if the Fe associated 
with the 44 Ti, it must be expanding in a direction close to the plane of the sky. The velocity 
dispersion indicates that the 44 Ti presumably occupies a volume significantly smaller than 
that outlined by the “pure” Fe, with a half angle subtended at the explosion center of order 
sin” 1 (430/4000) ~ 6 degrees. 


4.3. The Neutron Star Kick 

The neutron star is known to be recoiling from the explosion center inferred from the 
optical knot motions with a velocity in the plane of the sky of 330 km s” 1 (Thorstensen et al. 
2001). The direction of recoil is also nearly perpendicular from the NE-SW axis established 
by the jet and counter-jet, which presumably delineate an axis of rotation in the progenitor. 
This makes a kick arising from the action of a standing accretion shock instability appealing. 
It is known, at least in the limit of slow rotation, where distortion of the shock front from 
spherical symmetry can be neglected, that the instability grows fastest in the prograde 
spiral mode (e.g. Blondin & Mezzacappa 2007; Yamazaki & Foglizzo 2008; Iwakami et al. 
2009), i.e. , perpendicular to the symmetry axis. The pulsar kick imparted by the second 
supernova that formed the double pulsar system PSR B1913+16 is also inferred to have 
been largely perpendicular to the progenitor spin axis (Wex et al. 2000). The magnitude of 
the observed recoil velocity is also well within the accessible range generated by numerical 
simulations of such instabilities, albeit these are so far for initially nonrotating supernova 
cores (Wongwathanarat et al. 2010; Nordhaus et al. 2010). At higher rotation rates, it 
appears that the instability is suppressed (Burrows et al. 2007), though the deformation of 
the accretion shock by rotation may seed an axisymmetric quadrupole oscillation (Iwakami 
et al. 2009). 

In the case that the kick involves instability at the accretion shock, following the sus- 
picion advanced by Wongwathanarat et al. (2010), we would expect the Fe ejecta to recoil 
predominantly in the opposite direction to the neutron star. In Gas A, this appears not to be 
the case. Most of the Fe, and especially most of the “pure” Fe, is recoiling towards the east, 
within 90 degrees of the neutron star motion. However, the remnant as a whole, whose mass 
is dominated by O, is actually recoiling to the north, close to 150 degrees from the neutron 
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star motion, with a plane of the sky velocity of about 700 km s~ 4 , which perhaps suggests 
a hydrodynamic origin for the kick. Wongwathanarat et al. (2010) discuss a 15 M 0 explo- 
sion, whereas the Cas A progenitor was considerably less massive upon explosion such that 
reverse shocks and associated nucleosynthesis during the explosion are less likely to occur. 
Our estimate of the remnant recoil assumes homologous expansion, and the approximations 
made in assessing the ejecta mass have been detailed in sections 3 and 4.1. 


5. Conclusions 

We consider our efforts here to be a first attempt at a comprehensive view of the X-ray 
emitting ejecta in Cas A on few arcsecond angular scales. We have carried out at considerable 
effort the most detailed data analysis that was practicable for us, but there is ample room for 
improvement. The most obvious shortcoming of this work is the need for a more sophisticated 
definition of spectral extraction regions, and hand-in-hand with that, yet more extraction 
regions. In previous work, we had individually chosen regions to correspond to identifiable 
knotty or filamentary features. Here we were obliged for the sake of (relative) expediency 
to instead use regions that, while adjusted crudely for overall surface brightness, have been 
defined without reference to coherent structures in the surface brightness. Full use of the 
Chandra angular resolution in fact requires a more sophisticated definition of the extraction 
regions. The main advantage to be obtained, however, is that this may allow the use of 
simpler NEI models with a single ionization age, rather than parallel shock models with a 
linear distribution of ionization ages, and thus a much clearer connection with analytical 
hydrodynamical models (see for example, Laming & Hwang 2003; Hwang & Laming 2003). 
It may also eliminate the requirement for the filling factor correction used here. Finally, the 
analysis would also benefit from a more sophisticated treatment of the spectral background. 

We have arrived at a favored ejecta mass of about 3 M 0 based both on the observed SNR 
dynamics, and on an accounting of the X-ray emission and element abundances throughout 
the remnant. The increased signal/noise afforded by the long exposure time, and by the 
relatively large size of extraction regions has allowed us to identify regions of the remnant 
where the Fe ejecta is to a large extent “pure”, in that the mass ratio Fe/Si > 10. We seek to 
interpret this ejecta component as alpha-rich freeze out ashes, where the recently discovered 
(Iyudin et al. 1994; Vink et al. 2001; Renaud et al. 2006) 44 Ti and its decay products likely 
reside. The surprise here is that essentially all the Fe, and presumably also the 44 Ti and 
its decay products, are not located in the center of the remnant but have been ejected well 
out into the overlying ejecta, probably by the action of hydrodynamic instabilities during 
the explosion. In this respect, SN 1987A offers a stark contrast, with ejecta containing 44 Ti 
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located at the center with a characteristic expansion speed of rsj 3000 km s 1 (Kjaer et ah 
2010). Gl.9+0.3 may represent an intermediate case (Borkowski et al. 2010). 

Finally, Cas A has a well documented displacement of the neutron star formed in the 
explosion (Thorstensen et al. 2001), which is close to perpendicular to the NE-SW “jet”- 
axis. While contrary to expectations (e.g. Wongwathanarat et al. 2010), the Fe in Cas A 
does not recoil in the opposite direction to the neutron star, as far as we can tell, but the 
remnant as a whole does. We take this as reinforcing the notion that neutron star kicks in 
core-collapse supernova explosions may derive from purely hydrodynamic mechanisms, as 
recently suggested (e.g. Nordhaus et al. 2010). Together with motion and location of the Fe 
ejecta, and the likelihood that what is observed as a NE-SW “jet” really is an artefact of 
an asymmetric explosion, we speculate that an instability of the standing accretion shock 
(SASI) in a rotating progenitor is the most promising model for understanding the explosion 
that formed Cas A. 
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server Program and NNH04ZSS001N to the Long Term Space Astrophysics Program. JML 
was also supported by basic research funds of the Office of Naval Research. 


A. Revisions to Appendix A in Laming &; Hwang (2003) 

We collect here a few minor corrections and modifications to the text in Appendix A of 
Laming & Hwang (2003). Equations Al and A2 should read 

to = 473.6Mf/ 6 £ 5 - 1/2 p- 1/3 yr, (Al) 

xo = 3.43M e 1 /V 1/3 pc. (A2) 

The expressions for the lead factor and the pressure ratio between forward and reverse 
shocks following equation A5 are now more accurately given by 

8 0.4 / \ \ 

/ed — 1-0 H — - + (A3) 

n z 4 — s 

0ED = (0.65 - exp (— n/4)) \A - s/8. (A4) 

Finally, as pointed out to us by Elisabetta Micelotta, some of the discussion following 
equation A8 in Laming & Hwang (2003) requires clarification. Equation A5 can be derived 
from the envelope form of equation A8, i.e. from equations A7 and A8 in Truclove & McKee 



(1999). This procedure requires v e j = \J (10/3) Jn — 5) / Jn — 3) going from equation 74 
to 75 in Truelove & McKee (1999). In their equation 31, Trnelove & McKee (1999) 
give v core = \J (10/3) Jn — 5) / (n — 3), which implies w core = v core /v e j = 1. However their 
derivation of equation 31 assumes w core — > 0 in their derivation. A “cleaner” route to 
equation 75 of Truelove & McKee (1999) is to take their values of v core and the ejecta mass 
fraction in the core derived under the condition w core — > 0, and substitute into equation 3 of 
Chevalier (1982). The condition R b = /ED'iWeAore still gives the time t core when the reverse 
shock hits the ejecta core. The only other change connected with this discussion is that 
equation All for t conn in Laming & Hwang (2003) should depend on v core and not v e j. 

In actual fact, the difference between the approximations w core = 0 and w core = 1 is 
not all that large; v core = yj (10/3) (n — 5) / (n — 3) — > \/2 in these limits, which coincide for 
n = 8. The core fraction of the ejecta mass varies between (n — 3) jn and 1, and only v e j 
which goes from oo to v core between these two limits is sensitive. 
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Table 1: Model Parameters 


M ej /M & 

£51 

Vb 

Rb 

v r 

R r 

pr) 

F^bub 

^ exp 

M ej , x /M & 

m 

age 

^ core 

2.0 

1.7 

5376 

2.6 

2611 

1.65 

10 

0.26 

4900 

1.84 

0.694 

328.7 

10086 

2.5 

2.0 

5330 

2.6 

2509 

1.67 

12 

0.28 

4924 

2.28 

0.697 

332.6 

9785 

3.0 

2.35 

5357 

2.6 

2465 

1.70 

14 

0.30 

5012 

2.71 

0.699 

331.6 

9683 

3.5 

2.67 

5353 

2.6 

2416 

1.73 

16 

0.33 

5099 

3.13 

0.699 

331.9 

9555 

4.0 

3.0 

5357 

2.6 

2384 

1.75 

18 

0.35 

5163 

3.55 

0.699 

331.8 

9474 


Note. - Table columns give in order: M e j ejecta mass (M 0 ), E 5 1 explosion energy (10 51 ergs), Vb forward 
shock velocity (km s _1 ), Rb forward shock radius (pc), v r reverse shock velocity (km s -1 ), R r reverse shock 
radius (pc), pr 2 density times forward shock radius at the current position of the forward shock in H atoms 
(or equivalent mass)x pc 2 , Rbub circumstellar bubble radius (pc), v exp ejecta expansion velocity at reverse 
shock (km s -1 ), M e j >x mass of X-ray emitting ejecta (M 0 ), m expansion parameter, SNR age (years), and 
v core expansion velocity of ejecta core-envelope boundary (km s _1 ). 
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Table 2: Element Masses, Velocities 


element 

M e i t 2.5" / Mq 

M e i£"/M Q 

v x 

Vy 

M e i/ Mo 

M e l/ Mq\ solar 

O 

1.80 

2.55 

-176 

706 

1.00 

1.00 

Ne 

0.027 

0.038 

1232 

486 

0.015 

0.22 

Mg 

0.0070 

0.0099 

247 

323 

0.0039 

0.12 

Si 

0.038 

0.054 

-672 

677 

0.021 

0.12 

S 

0.020 

0.028 

-685 

787 

0.011 

0.053 

Ar 

0.010 

0.015 

-656 

920 

0.0056 

0.013 

F e Si 

0.073 

0.10 

-486 

420 

0.041 

0.23 

Fe a 

0.029 

0.041 

-1900 

-104 

0.016 


M x 

2.00 

2.84 

-210 

680 



M to t 

2.18 

3.14 

-210 

680 




Note. — Solar abundance ratios from Grevesse et al. (2010). Inferred masses scale as the square root of 
the filling factor. The Fe masses are taken for a probability threshold 0.01 as discussed in the text and given 
in Table 3. 












Table 3: Fe^, Fe Q Masses 


F e Si, 2 . 5 " 

Fe Qi 2.5" 

Fesi j5 » 

■P'^O',5" 

probability 

0.061 

0.095 

0.086 

0.13 

o 

o 

0.070 

0.049 

0.098 

0.070 

o 

1 

0.073 

0.029 

0.10 

0.041 

10~ 2 

0.075 

0.019 

0.11 

0.027 

10~ 3 

0.076 

0.012 

0.11 

0.017 

10~ 4 

0.076 

0.0092 

0.11 

0.013 

10~ 5 

0.077 

0.0063 

0.11 

0.0089 

10" 6 

0.077 

0.0050 

0.11 

0.0071 

o 

l 

-4 

0.077 

0.0043 

0.11 

0.0060 

o 

1 

00 


Note. — Fe masses from incomplete Si-burning (Fegj) and from complete Si-burning and/or a-rich freeze 
out (Fe a ) in M 0 . The final column gives the probability that the second Fe component is not justified by 
the reduced % 2 resulting from its inclusion in the fit. Demanding a lower probability reduces Fe a . 


